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ABSTRACT 

The disc of galaxies is made of the superposition of a thin and a thick disc. Star formation is hosted in the thin discs and contributes 
to their growth. Thick discs are formed of old stars. The formation mechanisms of thick discs are under discussion. Thick discs might 
have formed either at high redshift on a short timescale or might have been built slowly over a Hubble-Lemaitre time. They may have 
an internal or an external origin. Here we adopt a galactic archaeology approach to study the thick disc of ESO 533-4, i.e. we study 
the kinematics and the stellar populations of this galaxy in detail. ESO 533-4 is a Southern, nearby, and almost bulgeless galaxy. 

We present the first ever Integral Field Unit spectroscopy of an edge-on galaxy with enough depth and quality to study the thick disc. 
We exposed ESO 533-4 with the blue grism of the VIMOS instrument of the VET for 6.5 hours. The field of view covered an axial 
extent from ~ 0.1 r 2 s to ~ 0.7 r 2 s, where r 2 s is the 25 mag arcsec^^ isophotal radius. This corresponds to the range from ~ 1 kpc to 
~ 7 kpc. We used pPXF and the MILES library to obtain velocity and stellar population maps. We compared our kinematic data with 
simple GADGET-2 models. 

The apparent rotational lag of the thick disc of ESO 533-4 is compatible with that expected from the combinations of two etfects: 
ditferential asymmetric drift and the projection effects arising from studying a disc a few degrees (2 - 3°) away from edge-on. Thus, 
ESO 533-4 contains little or no counter-rotating material. This is compatible with three formation scenarios: the secular dynamical 
heating of an initially thin disc, the formation of the thick disc at high redshift in an early turbulent disc phase, and the creation of 
a thick disc in a major merger event. If this last mechanism occurred in all galaxies, it would cause retrograde thick discs in half of 
them. These retrograde discs have not been observed in the five massive disc galaxies (circular velocity Vc > 120kms“*) for which 
the kinematics of the thick disc is known. The stellar populations map indicates that the populations of the thin and the thick discs 
of ESO 533-4 are possibly separated in the Age - log (Z/Zq) plane. This would imply that thin and thick discs are formed of two 
distinct stellar populations. The stellar population results are not fully conclusive because of the high dust extinction in ESO 533-4 
and because recovering stellar populations is a difficult inverse problem. Having said that, the stellar population results do not favour 
a secular evolution origin for the thick disc. Hence, we suggest that the thick disc of ESO 533-4 formed in a relatively short event. 

Key words, galaxies: individual (ESO 533-4) - galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure - galaxies: 
evolution - galaxies: formation 


1. Introduction: Thick discs and galactic 
archaeology 

We kno w with absolute certainty that our uni verse contains 
galaxies (ICurtisl|1917l; IOpik|[l9^ iHubblel [19251) . Another fact 


established with almost the same degree of certitu de is that 


the initial state of the universe was a hot Big Bang (iLemaitre 

1931 

:lGamowll946l:lAlDher & Hermanll949l:lPenzias & Wilson 

1965 

). It follows that the products of the primordial nude- 


osynthesis somehow assembled and evolved into the galaxies 
that we observe roughly a Hubble-Lemaitre time after creation. 
The current c osmological paradigm is La mbda cold dark mat¬ 
ter (ACDM; iPlanck Collaboration! l2015l for a recent param¬ 
eter determination), where the energy budget of the universe 
is dominated by elusive dark matter and a positive cosmologi¬ 
cal constant. Within that paradigm, the evolution of galaxies is 
driven by a combination of interactions with the environment 


* Based on observations made at the European Southern Observatory 
using the Very Large Telescope under programme 09LB-0228(A). 


(nurture; e.g. Toomrel 1977]) and internal evolution (nature; e.g. 
iKormendv & Kennicutd 120041 lAthanassoulal 120131) . However, 
the details and the specific weight of each of these two families 
of processes in the galaxy construction and subsequent evolution 
still elude us. 

One of the most mysterious remnants of galaxy formation 
and evolution are thick discs. Seen as roughly exponential ver¬ 
tical excesses of light in edge-on galaxies a few thin disc scale 
heights above the galaxy mid-planes, their origin is still a mat- 
ter of debate more than three decennia after their discovery 
(iBursteinI Il979t iTsikoudil Il979h . Nowadays, all disc galaxies 
where thick discs have been sear ched for have been found to 
host one or even two of them ( [V^chim_& Dalcantc^ l2006t 
IComeron et aLll^l Idlbl but see IStreich et al.l (12015 )). The fol¬ 
lowing scenarios have been proposed to explain their formation: 


- External fast processes: the thick disc is the consequence of 
the merger of two or more gas-rich ga laxies during the initial 
assembly process (iBrook et al.ll2004l) . 
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Table 1. Summary of the signatures predicted in the thick disc formation scenarios 


Scenario 

External/ 

Internal 

Fast/ 

Secular 

Counter-rotating material 
in the thick disc 

Chemical and/or age discontinuity 
between the thin and thick discs 

Merger of gas-rich galaxies 

External 

Fast 

Possibly" 

Yes 

Turbulent early disc 

Internal 

Fast 

No* 

Yes 

Accretion of stars from satellites 

External 

Secular 

Possibly" 

Yes 

Heating by satellites 

External 

Secular 

No 

No 

Internal heating 

Internal 

Secular 

No 

No 


_ Those proces ses would cause some amount of counter-rotating material provided that some of the mergers are retrograde. In the 

iBoumaud et al.l d2014li simulations, some counter-rotating stars are created in this process. The fraction of retrograde stars decreases as a function 
of the mass of the galaxy. A high-z galaxy with a 1.4 x 10*® AIq haryonic mass, similar to the mass of ESO 533-4’s thick disc, would have 2.2% 
of counter-rotating stars (F. Bournaud, private communication). This small fraction of counter-rotating stars would he undetectable with the 

techniques used in this paper. 


Internal fast processes: the thick disc is born thick be¬ 
cause of the turbulent and clumpy nature of the first 
discs (Mmegreen & Elmegreenll20o5 iBournaud et al.ll^OOl: 


IComerdn et alj|20l4. 

External secnlar processes: the thick disc is made of stars 
stripped from infalling satellites dAbadi et al .1120031) and/or 
by the disc dynamical heating c aused by these satellites 
(lOuinn et al.lll993l : lOu et al.ll201 Ih . We call these processes 
secular because minor mergers are much slower than ma¬ 
jor mergers as a result of weaker dynamical friction. Also, 
according to ACDM, seyeral of these eyents can occur in 
a Hubble-Lemaitre time (e.g. seyeral tens of minor merge rs 
causing tidal features are expected in Llohnston et al.ll2008h . 
Internal secular processes: the thick disc is caused by dy- 
namical heating due to disc oyerdensit ies (IVillumsenlll985h 
and/or t he radial migration o f stars dSchonrich & Binn^ 
l2009allbL iLoebman et alJl 20 Tlh . Radial migration as a yiable 
mecha nis m to build thick discs is c ontested by lMinchey et alJ 
d2Q12h and IVera-Ciro et aP d2014l) . 


It is likely that a combination of two or more of these sce¬ 
narios contributed to the thick disc formation. Actually, some 
mechanisms, such as secular internal heatin g, must always oc¬ 
cur to some extent. In IComeron et alJ d2012l) . we suggested that 
thick discs in low- and high-mass galaxies might haye differ¬ 
ent dominant formation mechanisms. This is because there is a 
clear bimodality in the ratio of thick to thin disc masses as a 
function of the total galaxy mass; high-mass galaxies haye the 
same ratio irrespectiye of the galaxy total mass, whereas for 
low-mass galaxies the ratio decreases as the total galaxy mass 
increases. The fact that most high- mass galaxies haye rou ghly 
similar thick to thin disc mass ratios dComeron et alJl2()T4l) sug¬ 
gests the possibility that the same combination of mechanisms 
built them. The diyiding line between low- and high-mass galax¬ 
ies is_Jound_jt_a_£ircul^velocity_vy«i 120 kms“* (suggested 
by lYoachim & Dalcantonll2006ll2008bl) . which corresponds to a 
haryonic mass of Al ~ 10*° AIq when using the Tully-Fisher 
relation in IZaritsky et al.l d2014t) . 

A particular model dSchonrich & Binneyl l2009allbl : 
iRoskar et alJ 1201 iMinchey et al. 2015h . which proposes a 
thick disc created through a combination of internal and external 
secular heating mechanisms, is releyant in later sections of 
the paper. In this model, each of the mono-age populations 
of a galaxy is distributed in a flared disc. The younger stellar 
populations haye a longer scale length than the older stellar 
populations. The superposition of a series of flares at different 
radii produce what has been identified in photometric decompo¬ 
sitions as a thick disc. The result of this model is a galaxy with 


no apparent flare eyen if the indiyidual mono-age populations 
do flare. 

As a result of cosmological dimming and resolution prob¬ 
lems, it is often impractical to tackle galaxy eyolution issues by 
obserying highly redshifted objects. An alternatiye is so-called 
galactic archaeology, which consists of studying the structures of 
redshift z ~ 0 galaxies in great detail to infer their origin. Thick 
discs are ideal targets for galactic archaeology because differ¬ 
ent formation scenarios predict distinct kinematical and chemi¬ 
cal signatures, as summarized in Table[T] 

So far, attempts to do galactic archaeology in external 
edge-on galaxies to study thick discs haye been mostly based 
on photometric decompositi ons (lYoachim & DalcantonI l2006t 
Comeron_et_alJj20]Jj , l2012h. The exception is the work by 
Yoachim & Dalcanton (l2008bB) . in which single-slit spec¬ 
troscopy was used to study thin and thick discs. They found that 
thick discs are made of old stars and that thick discs of some 
low-mass galaxies contain counter-rotating stars. None of the 
three massiye galaxies in their sample had measurable amounts 
of retrograde stars. The problem that arises from single-slit ob- 
seryations is that one has to rely on photometric decompositions 
to And the range of heights where the thick disc dominates be¬ 
fore placing the slit. Ideally, we would like to be able to obserye 
galaxies in a wide range of heights, so no assumptions about the 
extent of the thick disc haye to be made a priori. 


2. ESO 533-4 

We present the first integral held spectroscopy study of 
the thin/thick structure of an edge-on disc. We used the 
Visible M ulti-Object Spec t rograp h of the Very Large Telescope 
(VIMOS; iLe Fevre et alJ 120031) . The selected target was 
ESO 533-4. This galax y has been part of the samp le of our re¬ 
cent thick disc studies dComeron et alJl^l2L l2014l) . The target 
was chosen because of: 

- its small central mass concentration (CMC), and hence, pre¬ 
sumably little effect of a central spheroid on the thin and 
thick disc kinematics; 

- its relatiyely small angular size, so a large fraction of the disc 
could be coyered by the VIMOS held of yiew (FOV); and 

- a brighter than ayerage thick disc, which starts to domin ate 
at m 3 6 ;jm(AB) 22 mag arcsec”^ dComeron et alJl201^ so 
that thick disc spectra can be obtained with a few hours of 
exposure. 

A summary of the properties of ESO 533-4 is presented in 
Table |2] There is a disagreement about the stage of the galaxy 
in the literature. ESO 533-4 has a small CMC, which makes it 
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Table 2. Parameters of ESO 533-4 


Parameter 

Value 

Source 

RA (J2000.0) 

22‘’14"'03’;02 

1 

Dec (J2000.0) 

-26°56'15'.'8 

1 

4/3.6/im(AB) 

-20.52 

2 

Type 

Sc? sp 

3 


S0° sp/E(d)8 

4 

Vsys (helicoentric) 

2596±3kms-‘ 

5 

D 

39.8 Mpc 

6 

1" physical size 

193 pc 

6 

Vc 

147.5 ±2.2 km s-' 

7 

r25 

5476 

7 

Luminosity profile 

Type 2 break 

8 

M, 

(680 ± 20) X 10’ Mq 

9 

Mcmc 

(19±2)x 10’ Mq 

9 

Ml 

(1500±300)x IO’Mq 

10 

Mt 

(1100±200)x 10 ’Mq 

10 


Refer ences, tli [ Skrutskie et all 1200^ : 12 1 iMuno z-Mateo s et all 
( l2015h : (3) Ide Vaucouleurs et all l li99lli : (4) iButaetalJ j2015l) ; (5) 
NED; (6) from NED average of r edshift -independent d is tances ; 
(7) HyperLeda dMakarov et all l2014h ; (8) IConieron et all ( 120121) ; 
(9) atomic gas disc mass (Mg) and CMC mass (Mcmc) from 
IComeron et all ( 120141) ; (10) thin disc mass (Mt) and thick disc mass 
(Mt) from this paper. 


appear late type (Sc: lde Vaucouleurs et aLlfToOlh . but its smooth 
appeara nce at 3.6 um ha s recently caused it to be classified as 
an S0° (iButa et al.ll2015l) . ESO 533-4 has a baryonic mass M ~ 
3 X10A4 q , and hence it is a high-mass galaxy. A 3.6 /rm image 
from the Spitzer Survey of Stellar Structure in Galaxies (S'^G; 
ISheth et al.ll2010l) is presented in Eig. [1] 


3. The thin/thick disc decomposition of ESO 533-4 


In IComeron et al.l (1201 2h . we made thin/thick disc decomposi¬ 
tions of ESO 533-4 for the symmetrized 3.6/rm + 4.5/rm lu¬ 
minosity profiles in four bins perpendicular to the galaxy mid¬ 
plane. The discs were assumed to be in hydrostatic equilib¬ 
rium following the formalism bv iNaravan & .lod (l2002h . We ac¬ 
counted for the gravitational effect of a thin gas disc with 20% of 
the surface mass density of the thin disc. We made several sim¬ 
plifying assumptions, such as similar scale lengths for the thin 
and thick discs, that the disc is not too submaximal in its inner 
parts, vertical isothermality, and roughly constant scale heights 
with c hanging radii. 

In IComeron et al.l (1201 2h we convolved our solutions with a 
Gaussian kernel to account for the point spread function (PSE) 
of the S^G images. The kernel had a 2"2 full width at half max¬ 
imum (EWHM). However, the S^G PSE is now known to have 
faint but extended wings not accounted for by a Gaus sian. Here 
we re do the fits using the S'^G symmetrized PSE (ISalo et al.l 
l2015l) . Here we produce our fit using only the 3.6yum band be¬ 
cause we only have a detailed PSE model for that band. 

We use the sky va lues from the S^G Pipeline 3 
dMunoz-Mateos et al .1120151) instead of our own sky determina¬ 
tions. To do the fit, we use idl’s c urvefit instead of our own 
m inimization rou t ine (d escribed in IComeron et al.ll201 lah . As 
in IComeron et al.l (l2Q12l) . we use a ratio of mass-to-light ratios 
of the thick and thin discs Tx/Tt = 1.2. The Tj/Tt = 1.2 
value comes from studying the spectral energy distributions 
res ulting from the Milky Way st ar formation history models 
by iNykvtvu k & Mis henina] (l2006l) for the thin and the thick 
discs ( Comeron et al.l 201 lah . The maximum height of the fit 



Fig. 1. 3.6 yum image of ESO 533-4 from the S"^G. The red cross 
indicates the centre of the galaxy and the region with a green 
border is that we studied. North is up and east is left. The white 
solid lines perpendicular to the galaxy mid-plane show the four 
axial bins used for the luminosity profile fits in Eig.|2 The dash- 
dotted lines indicate the heights over which the profiles were fit. 
The solid lines parallel to the galaxy mid-plane correspond to 
the height above which 90% of the light comes from the thick 
disc. The numbers are used to associate the bins with the fits in 
Fig.m 


was selected as that beyond which it is not possible to make 
a two disc component fit (the mean squared difference be¬ 
tween the observed profile and the fit becomes larger than 
0.01 (mag arcsec^^)^ at that height). What limits the height of the 
fitting region is the presence of scattered light of point sources in 
the field and background gradient issues. The extent of the fitting 
range covers over four thick disc scale heights. 

The new thin/thi ck disc fits are very similar to those in 
IComeron et al.l(l2012l) . which indicates that the errors introduced 
by approximating the PSE with a Gaussian function were small. 
The symmetrized PSE used here covers 15" in radius. The thick 
disc dominates the luminosity profiles down to z 5" (~ 1 kpc), 
which is a height safely smaller than the radius of the modelled 
PSE. If we consider our fitted thin disc, we find that the extra 
light added to |z| > 5" is about 4% of the observed emission. 
Therefore, we can discard the possibility that the thick disc is a n 
artefact caused by scattered light as suggested bv lSandinl(l2015l) . 


The new fits indicate a slightly higher ratio of surface den- 
sities between the th ick and thin disc (Sj/St) compared to 
IComeron et al.l (1201 2l) . which slight ly increases the mass ratio 
of the two discs: Mx/Alt = 0.65 in IComeron et al.l (120121) an d 
My!M l - 0.77 here, following Eq. 5 in Comeron et al. ( 201 lal) . 
We could use a gas disc surface density larger than 20% of that 
of the thin disc to reflect the fact that the gas mass is ~ 45% that 
of the thin disc (Table |2](. However, this would not take into ac¬ 
count that usually gas discs are more extended than stellar discs. 
Furthermore, when we refitted the luminosity profiles with dou- 
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Fig. 2. Luminosity profiles of ESO 533-4 in the four bins perpendicular to the mid-plane that are indicated in Fig.[T](large dots). The 
fits are indicated with the dotted lines, and the contribution of the thin and thick discs are indicated with dashed lines. The vertical 
solid line indicates the height, Zs, above which 90% of the light comes from the thick disc according to our fit. The dot-dashed grey 
vertical line indicates the maximum height, z\, that we use for the fit. pro/Pto stands for the thick to thin disc mass density ratio at the 
mid-plane, CT'i/crt for their velocity dispersion ratio, zj and Zt for the thick and thin disc vertical scale heights, Sj/St for the thick 
to thin disc mass ratio, pi for the lowest fitted surface brightness level, and Am for the range in magnitudes over which the fit was 
done. 


ble the gas surface density, we obtained the same Afr/Aft within 

1 %. 


We find that 90% of the light coming from the region |z| > 5" 
(z > 970pc) is emitted by the thick disc. In what follows, we 
consider that region as that where the thick disc dominates the 
surface brightness. From our f its, we obta i n that 39% of the light 
comes from the thick disc. In Salo et all (l2015h . this fraction is 
33% when using Galfit3.0 (iPeng et al.ll2010h for the decompo¬ 
sition (Galfit3.0, though performing a simultaneous 2D fit over 
the whole galaxy, oversimplifies the luminosity profiles perpen¬ 
dicular to the mid-planes by assuming that the discs have verti¬ 
cal sech^(z/zo) profiles). The Galfit3.0 thick disc scale height is 


slightly underestimated (zx = 3'.'4 for Galfit3.0 versus zt - 4'.'4 
in our determination when averaging over the four fits in Fig.|2]i. 
Galfit3.0 finds that the height above which 90% of the light 
comes from the thick disc is z = 4'.'6. 


4. Observations 

The spectroscopic observations were made on the nights 6-7 
and 7 - 8 of August 2013. We used the integral field unit mode 
of VIMOS with the High Resolution Blue grism and a spacial 
resolution of Ol'hh x (y.'66 per fibre. The FOV of this mode is 
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27" X 27". The wavelength coverage of this conhguration is 
3700 - 5350 A and the dispersion is 0.71 Apixer*. 

The total in-target exposure was 6.5 hours (2 hours in the 
hrst night and 4.5 hours during the second night). The exposure 
was divided into 13 observing blocks (OBs) with a few arcsec- 
ond dithering between them. Each OB exposed ESO 533-4 for 
30 minutes. We also took 4 minute exposures of the sky, at a 2' 
distance in the direction perpendicular to the mid-plane. During 
the first night, a sky exposure was taken for every two galaxy 
OBs. During the second night two sky exposures were taken 
for every three galaxy OBs. Both nights were photometric and 
dark. The seeing was large (> 1"), but this does not affect our 
science because in any case we need a large spacial binning to 
study the thick disc with a sufficiently large signal-to-noise ratio 
{S/N). Standard calibration exposures (flat fields, bias, calibra¬ 
tion lamps) and spectrophotometric standard stars were taken for 
the two nights. 


5. Data reduction 

5. 1. Basic reduction 

Eirst, we used idl’s la_cosmic routine to clean the cosmic rays. 
The spectra were reduced and flux calibrated by the VI MOS 
pipeline. The Reflex environment dEreudling et al.l l2013h was 
used to run the pipeline. VIMOS has four quadrants and the 
pipeline provides independent reductions for each quadrant. 
Therefore, the output of the initial reduction was 4 x 13 = 52 
sets of spectra corresponding to the galaxy. 

Eor each of the quadrants of the sky exposures, we created 
a single spectrum by averaging the spectra in each spaxel with 
a 3cr clipping. We calculated the sky emission for each of our 
galaxy exposures by associating them with the average of the 
sky spectra from the same quadrant taken within 64 minutes. 
This time was set so at least one sky exposure could be linked to 
each galaxy exposure, but in some cases there were up to four. 
Each of the quadrants has a remarkably uniform sky background. 
Quadrant number 2 has ~ 40% less signal than the others. 

The information in headers could not be used for spacial 
alignment because of a few arcsecond imprecisions in the point¬ 
ing. To circumvent this problem, we tried to align the frames by 
hnding point sources in images made by collapsing the spectra. 
However, the two point sources in the vicinity of the pointings 
were close to the edge of the field, and thus did not show up 
in some of the exposures as a result of the ditherings. Then, 
we built [Out] images by integrating in the spectral direction 
from A - 5047.5 A to 4 = 5054.6 A (this takes the redshift of 
ESO 533-4 into account). Unresolved star-forming regions stand 
out in the [Out] images and were used for manual alignment. 
Once the galaxy frames were spacially aligned, we summed the 
spectra for every spaxel using a 3cr clipping. We also summed 
for each spaxel the sky spectra associated with the galaxy frames 
that were not clipped away. Eor every wavelength in each spec¬ 
trum corresponding to each spaxel, we estimated the S/N by 
considering the Poisson noise and the readout noise, 

S/N = (1) 

yjs galaxy + 5 sky + 2; {RON,^^ 


S/N because it is very small (typically 5 hr'' pix'*, according 
to the VIMOS Pipeline User ManuaQ)- 

We masked the two point sources in the held, thought to be 
foreground stars or ESO 533-4 globular clusters, as well as sev¬ 
eral spaxels in the edges that were covered by few exposures. 
The two point sources in the held were used to align the lEU 
data with the S^'G data. The alignment precision is estimated to 
be of the order of 1". 


5.2. Kinematics 


5.2.1. Velocity maps 

Spacial binning is necessary to achieve a S /N good enough to 
study thick disc kinemati cs. The binning wa s done with the 
Voronoi binning code by ICappellari & Coninl (l2003h with the 
condition of S /N ~ 25 for spectral pixels between A - 4500 A 
and A - 5000 A. 

To hnd the velocity held of ESO 533-4 we htted the spa¬ 
cially binned spectra with the MIL ES stellar population syn¬ 
thesis models (IVazdekis et^ 2010h using the penalized pixel- 
htting (pPXE) code by Cappellari & EmselletiJ (l2004lPl It is 
usually required to Gauss-convolve either the library or the ob¬ 
served spectra so their EWHMs match. However, we found that 
the EWHM in our spectra is AT « 2.2 A on the blue side and 
AT 2.8 A on the red side, which is very similar to MILES’ 
EWHM, AT = (2.51+0.07) A. We therefore decided not to make 
a convolution. Because of the relatively low S/N of the bins, 
we only htted two momenta. Emission lines were masked with 
a mask Av = 1000 km s'' in width. We used an eight degree 
additive Legendre polynomial to correct the continuum shape 
for dust reddening and possible photometric calibration prob¬ 
lems. The velocity map is shown in the right panel in Eig.[3] The 
height above which 90% of the light comes from the thick disc 
is |z| 5". Therefore, the uppermost and the lowermost rows of 

bins should be dominated by the thick disc light. 


5.2.2. Rotation curves 

We made rotation curves of ESO 533-4 at hve different heights, 
namely z = 0, z - ±3" (~ 0.6 kpc), and z = +8" (~ 1.5 kpc; 
Eigure|4]i. The z = ±8" cuts should be well within the region 
where the light emission comes mostly from the thick disc. The 
prohles take into account the movement of the Earth away from 
the galaxy (~ 19 km s ' at the time of the year at which the ob¬ 
servations were performed) and assume a heliocentric systemic 
velocity Vsys = 2596 km s ' (Table|2]i. The points in the rotation 
curves are calculated by making axial cuts in the data used for 
Fig .[3] The error bars in the rotation curves are estimated by pro¬ 
ducing 20 Monte-Carlo simulations (with the estimated noise for 
each pixel) when obtaining the kinematics from each spectrum. 

When extrapolated to R = 0, the rotation curves are close 
to V = 0, which conhrms that the spectra calibration process is 
robust. It is also reassuring that our z - 0 rotation curve matche s 
very well the Ha rotation curve from iMathewson et al.l (Il992h . 
once this has been offset by +10kms''. The correction seems 
necessary because taken as it is published, the velocity at R = 0 
is about -10km s '. If the offset is not made, the maximum in 
the Ha rotation curve is found at v ~ 160 km s ' on the west- 


where, 5 galaxy is the signal from the galaxy, 5 sky is the signal 
from the sky, and RONi are the read-out noises of each of the 
exposures that were summed to obtain the hnal spectrum of a 
given spaxel. We neglected the dark current contribution to the 


' ftp://ftp.eso.org/pub/dfs/pipelines/vimos/vimos-pipeline-manual- 
6.8.pdf 

^ Both pPXF and the Voronoi tesselation codes can be found at 
http://www-astro.physics.ox.ac.uk/~mxc/software/ 
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Fig. 3. Left panel; symmetrized, mid-infrared S^G luminosity profile and thin/thick disc decomposition of the 0.08 r 25 - 0.69 r 25 
(0.8 - 7.3 kpc) axial range of ESO 533-4 (the one covered by the IFU). The horizontal solid lines indicate the height, Zs above 
which 90% of the light comes from the thick disc. Middle panel: image of the observed held made by collapsing the spectra. Right: 
velocity held of ESO 533-4 after the subtraction of the recession velocity. The horizontal axis indicates the axial distance to the 
galaxy centre and the vertical axis indicates the distance from the mid-plane. 
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Fig. 4. Rotation curves of ESO 533-4 at different heights as in¬ 
dicated by the colours. The northern half of the disc corresponds 
to the bottor n half of the disc in Fig . |3] The black dots corre¬ 
spond to the iMathewson et al.l (119921) Ha rotation curve with a 
- 1-10 km s“' offset (see text). 

ern side of the galaxy and at v ~ 140km s“' on the eastern 
side, which would be strange in an undisturbed canonical rota¬ 
tion curve. We suspect an inaccurate calibration or heliocentric 
correction to be the source of the problems with the Ha data as 
published. 

5.3. Stellar populations 

To study the stellar populations, we manually created new spa- 
cial bins by combining those used in the velocity map. Before 
summing all the spectra from the low-S/N bins to create the 
high-5 /A^ binning, we shifted their wavelengths so they all corre¬ 
sponded to a 2500 km s ' recession velocity. We made the bins 
roughly rectangular, with the long side in the direction of the 


plane of the galaxy, so different heights could be studied. To 
achieve the required S/N for the thick disc, we summed the 
light from the two sides of the mid-plane. Manual binning is 
essential because the Voronoi tesselation tries to make “com¬ 
pact” or “round” bins (ICappellari & Conhil120031) . This second 
binning yielded spectra with a S /N ~ 80 in the A = 4500 A to 
A = 5000 A spectral range. 


To obtain the stellar population map, the spectra were again 
htted using pPXF. We used a grid of MILES stellar population 
templates ranging from 63Myr to 17.8 Gyr (50 age bins) and 
from log (z/Zq) = -1.71 to log (z/Zq^ = 0.22 (six metal- 
licity bins). Thus, the total template number was 300. A regu¬ 
larization factor of REGUL=200 was used to smooth the stellar 
populations. The continuum was corrected with multiplicative 
Legendre polynomia (for details on the regularization and the 
use of Legendre corrective polynomia, see the pPXF documen¬ 
tation). These corrections are essential because of the dust red¬ 
dening close to the mid-plane and because the blue side of the 
spectrum has little signal, so it is prone to photometric calibra¬ 
tion errors. We tested which polynomial order would work best 
for our observations. The thick disc ht, where dust reddening 
is likely to be negligible, is very stable with respect to changes 
of the order of the polynomial (which indicates little dust af¬ 
fecting the spectra and also a good data calibration that makes 
correcting polynomia unnecessary). The htted thin disc stellar 
populations are very dependent on the order for orders smaller 
than ~ 10. For larger orders the htted stellar populations change 
little with changes in order, which probably indicates that the 
continuum shape is properly corrected for differential dust ab¬ 
sorption and no further orders are required. The results presented 
here are htted with an order 14 polynomial, which is compara- 
ble to what others have used with VIMOS data (for example 
IChilingarian et al.ll2009l used 25). The stellar population map is 
presented in Fig. |5] The ht to the spectra in two of our tiles are 
shown as an example in Fig.|^ In Fig.|7] mean age and metallic- 
ity maps are shown. 
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Fig. 5. Stellar populations of ESO 533-4. The background image is the same as in the middle panel of Fig. [3] The red lines indicate 
the spacial bins used for the analysis. The horizontal solid yellow lines indicate the height above which 90% of the light comes from 
the thick disc. On top of each bin there is a stellar population plot with the horizontal axis corresponding to the ages in Gyr and the 
vertical axis corresponding to the metallicities in log (z/Zq^. The two thick disc sections (uppermost and lowermost tiles) were 
actually treated as a single bin. The colours in the plots for each bin indicate the mass fraction of a given stellar population. The 
bottom left corner of the figure shows an enlarged version of one of the stellar population distributions to improve the readability of 


the axes. 

6. Interpretation of the data 

6.1. Kinematics 

We find an asymmetry in the velocity map: within the thin disc 
at a given radius, the southern (upper) half of the galaxy seems 
to rotate faster than the northern (lower) half (Fig. [3 and the 
z - +3" or z = +0.6 kpc rotation curves in Fig. |4]i. Possible 
reasons for the asymmetry are the presence of a dust lane in 
a not perfectly edge-on disc and features that are not symmet¬ 
ric with respect to the mid-pl ane (e.g. some phases of a buck - 
ling bar such as those seen in iMartinez-Valnuesta et al]l2006h . 
The latter possibility, however, is unlikely because there is no 
sign of a strong bar in ES O 533-4 (boxy/peanut inner isophotes; 
ICombes & Sanderilll98lh . 


The rotation curves for the northern half of the galaxy (below 
the mid-plane in Fig. |3]l decrease monotonically in amplitude as 
z increases. For the southern half, this in only true for R > 25" 
(a striking feature of the rotation curves is that for the southern 
half of the galaxy, we find the velocities at z = 3" to be larger 
than those at z = 0 in the 12" -25" or 2.3 -4.8kpc axial extent). 
The difference in velocities between the mid-plane and |z| = 8" 
(z = ± 1.5 kpc) is 10-30 km s“' at the radius where the maximum 
circular velocity is reached. The apparent rotational lag of the 
thick disc must at least partly be a consequence of a stronger 
asymmetric drift caused by the larger velocity dispersion of the 
thick disc. Part of the lag might be caused by the presence of 
counter-rotating stars. Finally, part of the apparent lag might be 
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Fig. 7. Mean stellar age (left) and metallicity (right) maps. The white lines indicate the tesselation. The horizontal solid green lines 
indicate the height above which 90% of the light comes from the thick disc. 


because we may be looking at a disc that is not perfectly edge- 
on. 


To check for the causes of the asymmetries in the velocity 
map and for those of the l ag of the thick disc, we ran a set of 
models using GADGET-2 (ISt)ringelll2005h . The initial models 
were created with the GalactICS software dKuiiken & Dubinskil 
I 1995 I 1 . They included a thin disc and a thick disc, which had sim¬ 
ilar scale lengths. The scale height of the thick disc was ~ 4.9 
times larger than that of the thin disc and the mass ratio was 
1.1 /1.5 in agreement with the values in Table |2] and the fits in 
Fig. 12 Each of the discs was composed of 1.03 x 10® particles 
with a Toomre parameter Q - 1.3. We also included a small 
spherical CMC, with a mass 50 times smaller than that of the 
sum of the two discs, formed of 2 x 10"^ particles. Finally, we 
embedded the whole system into a dark matter halo made of 
1.08 X 10® particles. The distribution of dark matter particles was 
set to make a roughly flat rotation curve at large radii. At a ra¬ 
dius of two scale lengths the disc contributed to about 40% of 
the total radial force. The baryonic mass of our model scaled to 
the observed axial and vertical (Figure |9]) prohles and the rota¬ 
tional velocity amplitude is 2.9 x 10'® Mq. This is close to the 
actual mass of ESO 533-4. Models with different dark matter 
fractions could be created, but then having a good simultaneous 
vertical and axial match between the models and observations 
would not have been possible. The model snapshots shown here 
are for models that have been run for over ten rotations at all the 
radii covered by the VIMOS observations. This guarantees that 
the thin and thick discs have settled to equilibrium. 


To account for the effect of dust, we assumed it would follow 
the distribution of thin disc stars except for the fact that we con¬ 
sidered a scale height that was three times smaller. We arrived at 
the ratio of three times smaller because we assumed dust to be 
mixed with dense ga s. In the Milky Way , dense gas has a scale 
height of ~ 100 pc dLanger e t alJ 201 4ll wherea s the thin disc 
scale height is ~ 300 pc (IGilmore & ReidllToS^ . We assumed 
the dust would cause a grey absorption, and we studied several 


dust densities. The results shown here correspond to an optical 
depth of T = 8 to the centre of the galaxy when looking at it from 
a perfectly edge-on view. 

We checked whether ESO 533-4 has some measurable devi¬ 
ation from an edge-on orientation by looking at its colour map 
(Fig. EJ. We found the redder sections of the thin disc to lie 
slightly below the mid-plane (~ 1"). This implies that the near 
side of the galaxy is on top, which corresponds to the southern 
side of the observed galaxy in Fig. [3] We find that a ~ 1" devia¬ 
tion of the dust lane corresponds to a galaxy inclination i 85° 
in our models. We simulated colour maps by comparing the lu¬ 
minosity profiles of models with and without dust. ESO 533-4’s 
inclination, however, is likely to be closer to 90° given the ex¬ 
tremely thin aspect of the disc and the lack of evidence for spiral 
features (Fig.[T]). If the dust distribution were less centrally con¬ 
centrated than in our models, deviations of only 2-3° away from 
edge-on would be enough to explain the vertical shift of the dust 
lane. There are clear examples of galaxies where the dust dis¬ 
tribution is more extended than that of t he stellar disc (e.g. th e 
Sombrero Galaxy or the system studied in iHolwerda et alj2009ll . 

To create synthetic velocity maps, we assigned the MILES 
spectrum of a 2.5 Gyr star with solar metallicity to thin disc par¬ 
ticles. For the thick disc, we used the spectrum of a 10 Gyr star 
with metallicity Z = 0.2 Zq. We created synthetic images of the 
galaxy by projecting them and convolving them with a Gaussian 
PSF with a FWHM of 1'.'2. We then used the Voronoi tesselation 
software to create a tesselation matching that obtained for the 
VIMOS data (i.e. the density of tiles in the thin disc is similar). 
We then obtained the velocity maps exactly as we did for the 
VIMOS data. 

Figure [10] shows the velocity maps for several models in 
an area approximately matching that of the VIMOS observa¬ 
tions. The top row in the figure corresponds to galaxies with 
no dust and the bottom row corresponds to galaxies with dust. 
Models with no dust are symmetric with respect to the mid- 
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Fig. 10. Top: observed velocity map of ESO 533-4. Middle and bottom rows: set of velocity maps from models designed to mimic 
the region shown in Fig. [3] The three columns correspond to: a) a perfectly edge-on galaxy with no counter-rotating material; b) a 
galaxy 3° away from edge-on with no counter-rotating material; and c) a galaxy 3° away from edge-on with 10% of counter-rotating 
material in the thick disc. The top row corresponds to models with no dust and the bottom row corresponds to models with some 
amount of dust (see text). The horizontal lines indicate the height above which ~ 90% of the light comes from the thick disc. 


plane. Models with dust, when not seen exactly edge-on, have 
some degree of asymmetry. 

The left column in Fig.fTOlcorresponds to models with a per¬ 
fectly edge-on view of the galaxy. The optical depth is the largest 
close to the mid-plane. Dust attenuates the light from the point 
which is closer to the galactic centre (that with the largest veloc¬ 
ity along the line of sight) more than that of stars closer to the 
observer. Thus, the dust attenuation causes an apparent slowing 
of the stars around z = 0. 

The middle and the right columns in Fig. [10] show that dust 
is able to create asymmetries when looking at galaxies that 
are not perfectly edge-on (inclination i - 87° in the example 
shown here). The larger attenuation on one side of the galaxy 
causes the asymmetry in the velocity maps: the maximum ve¬ 


locity at a given radius is not found at the mid-plane, but slightly 
above it. The deviations from symmetry in our models with 
i = 87° and dust are in qualitative agreement with the obser¬ 
vations (Figure O. 

To find the cause of the apparent lag in the thick disc of 
ESO 533-4 (asymmetric drift, galaxy inclination, and/or pres¬ 
ence of counter-rotating stars), we created rotation curves at dif¬ 
ferent heights for the models with the same approach we used for 
the VIMOS data (Figure [TT]). We find that at the radius where 
v(z = 0) ss 150 km s“', the difference in velocity between the 
mid-plane and z = +8" (z = +1.5 kpc) is ~ 20kms“* for a 
model with no counter-rotating material and ~ 30 km s ' for a 
galaxy whose thick disc has 10% of counter-rotating material. 
These numbers are similar to those observed in ESO 533-4 and 
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Spectrum from the region dominated by the thick disc light 



I (A) 


Spectrum from a region dominated by the thin disc light 



I {A) 


Fig. 6. Two examples of rest-frame spectra obtained from our 
S /N ~ 80 tesselation. The black line corresponds to the actual 
spectrum, the red line corresponds to the fit, and the blue verti¬ 
cal lines indicate the spectral windows that were masked because 
of the possible presence of emission lines. Residuals are shown 
with green symbols below the spectra. The top spectrum belongs 
to the thick disc dominated region and the bottom spectrum to 
the rightmost thin disc dominated tile in Fig. |5] The spectra are 
normalized to their median value. Several spectral lines are indi¬ 
cated with dashed lines. 
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Fig. 8. Left: VIMOS-3.6//m colour map. The VIMOS im¬ 
age was degraded to the S"^G resolution. Right: median 
VIMOS-3.6jum as a function of the height with an arbitrary 
zero point. The solid horizontal lines indicate the height above 
which ~ 90% of the light comes from the thick disc. The dashed 
horizontal line indicates the reddest height within the thin disc 
dominated region. 
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Fig. 9. Luminosity profiles of the observed galaxy (large black 
dots) and the model galaxies (grey lines) in the 0.08 r 25 -0.69 r 25 
axial range (dust extinction set to zero). Two models, one with 
i - 90° (dash), and one with i - 87° (dots), are shown. The zero 
point of the model profiles is arbitrary. 


indicate that the thick disc has little or no counter-rotating mate¬ 
rial. The effect of the asymmetric drift plus the disc orientation 
(slightly away from edge-on) are enough to explain the thick disc 
lag. 

In Fig. [TT] we see how, in the southern half, the z = 3" 
(z - +0.6 kpc) curve overlaps with that at z = 0 for some ra¬ 
dial extent. This is because dust absorption is larger at z = 0 
than at z = 3". By increasing dust absorption it is possible to 
obtain the z = 3" curve above the z = 0 curve for some axial 
extent just as in the observations. This effect is therefore most 
likely due to dust. 

Our models are not able to reproduce the rotation curves of 
the northern half of ESO 533-4 forR < 25" (R < 4.8 kpc), where 


the thick disc lags by 30 km s“' or more with respect to the mid¬ 
plane. The mismatch between the thick disc and the z = 0 ro¬ 
tation curves could be due to some fraction of counter-rotating 
stars in the inner parts of the galaxy. However, if that were the 
case it would be difficult to explain why the counter-rotating ma¬ 
terial shows up only on one side of the galaxy. 

6.2. Stellar populations 

We distinguish two distinct stellar populations in the subplots in 
Fig. |5] We find: 

- an old (~ 10 Gyr) metal-poor population and 
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Fig. 11. Rotation curves of models at different heights as indi¬ 
cated by the colours. The northern half of the disc corresponds 
to the bottom half of the disc in Fig. [TOl The top row rota¬ 
tion curves correspond to a model with i - 87°, dust and 0% 
of counter-rotating material, and the bottom row corresponds to 
similar model with 10% of counter-rotating material in the thick 
disc. 


- a relatively young metal-rich population. 

Also, traces of young metal-poor populations, either a fit artefact 
or associated with the outskirts of the thin disc, are seen in a few 
of the tiles (this is most evident in the tiles 2, 4, and 6). 

To quantify the fraction of mass associated with each of these 
stellar populations in each tile, we divided the diagrams in two 
regions, as seen in Fig. [12] The limit between the thin and the 
thick discs is set at log (Z/Zq) = -0.7 for ages larger than 3 Gyr. 
All populations younger than 3 Gyr are considered to belong 
to the thin disc. While a f = 3 Gyr stellar population is young 
compared to what one would expect for a thick disc population, 
the regularization introduced by pPXF causes stellar populations 
with peak ages at f ~ 10 Gyr to extend down to smaller ages. 

The first of the populations in the list is naturally associated 
with the thick disc, and the second is likely to be related to the 
thin disc. To further confirm this, we find that the old metal- 
poor population dominates the mass budget in the bin with the 
largest height, whereas the younger metal-rich population dom¬ 
inates close to the mid-plane. 

In many tiles, the stellar population fits show that the thin 
and thick disc populations have two distinct maxima in the 
Age - log (Z/Zq) plane. In these tiles, the thin and thick disc 
stars seem to belong to two distinct populations. In other tiles 



Age (Gyr) 


Fig. 12. One of the subplots in Fig.|5](the rightmost). The green 
lines indicate the divisions between the thick disc population 
(bottom right) and the thin disc (rest of the plot). 
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Fig. 13. Ratio between the masses of the thin and thick disc 
stellar populations (top panel) and its inverse (bottom panel), as 
a function of the height |z| of the Fig.|5]tile. The lines indicate the 
same ratios according to the fit to the light profiles (Figure|3]l. 


(particularly 4, 6, 11, and 12 in Fig.|5]l the two populations are 
not so clearly separated. The choice of the regularization param¬ 
eter introduces a certain arbitrarity in the interpretation of the 
age-metallicity plane. However, this is a necessary evil, since a 
certain amount of regularization is needed to smooth the other¬ 
wise even more arbitrary solutions of a difficult inverse prob¬ 
lem (see pPXF online documentation). While varying the reg- 


11 


































Comeron, S., et al.: Galatic archaeology of the ESO 533-4 thick disc 


ularization parameter between REGUL=100 and REGUL=250 
does not drastically change the results, higher values can wash 
out significantly the peaks. The chosen regularization value of 
REGUL=200 is similar to the value REGUL=250 used in the 
only source where, to our knowledge, a concrete value of this 
parameter has been explicitly stated (the example codes from 
the pPXE website). 

If we assume that the stars are divided into a thin and thick 
disc population (as indicated by many of the Age - log (Z/Zq) 
plots in Eig. |5]), we can test the assumptions made when fitting 
the luminosity profiles t o decompose them into th e thin and thic k 
disc components, as in IComerdn et al.l (l201 Itfl 1201 2i l2014ll . 
We calculated for each tile the mass that corresponds to the thin 
and thick discs according to the fit in the left panel of Eig. |3] 
To accomplish this, we weighted the spaxels accounting both for 
the vertical and the axial lumino sity profiles (axial luminosity 
profile from IComeron et ^l2012h . Because of the broadness of 
the peaks that correspond to the stellar populations of the thin 
and the thick discs (see Eig. [nj, the thick disc stellar popula¬ 
tion might extend to larger metallicities that the ones consid¬ 
ered here (limit set at log(Z/ZQ) = -0.7). Likewise, the thin 
disc stellar population could extend to smaller metallicties than 
the ones considered here. To account for this, we considered 
the possibilities that the thick disc population extended up to 
log(Z/ZQ) = -0.3 or only up to log(Z/ZQ) = -1.1 (for ages 
larger than 3 Gyr). We also considered the possibilities that the 
thin disc extended either down to log(Z/ZQ) = -1.1 or only 
down to log (Z/Zq) = -0.3 for ages larger than 3 Gyr. This is 
what has been used to calculate the error bars in Eig. [T3] The 
figure shows that the agreement between the results obtained 
from the stellar population maps and that from our thin/thick 
disc modelling is good, especially if we take the uncertainties 
and assumptions made in the modelling process into account. 


7. Summary of the results, discussion, and 
conclusions 

In this paper we have made an lEU study of the thin and the thick 
discs of ESO 533-4 with VIMOS. The main findings are: 

- The velocity field of ESO 533-4 is compatible with that of 
a galaxy inclined a few degrees from edge-on with a signifi¬ 
cant dust absorption in the equatorial plane. The near side of 
the galaxy is the southern one. 

- The apparent velocity lag of the thick disc can be explained 
by asymmetric drift and because of the projection effects 
arising from studying a disc a few degrees (2 - 3°) away 
from edge-on. Eew or no counter-rotating stars are needed to 
account for the asymmetric drift of the thick disc. 

- The thick disc is made of old (~ 10 Gyr) and metal-poor 
stars. The thin disc is made of younger metal-rich stars. This 
is manifested as a distribution with two distinct peaks in the 
Age - log (Z/Zq) plane for the majority of the regions of 
ESO 533-4 we studied. This might indicate that the stellar 
populations of the thin and thick discs are distinct. In several 
regions (at least 4 out of 20) there is no distinct peaks, but 
instead a continuum distribution of stellar populations from 
old metal-poor stars to younger metal-rich stars. 

- The relative fractions of mass in the thin and the thick discs 
at different heights, according to our stellar population anal¬ 
ysis, agree with those obtained from our thin/thick disc fits 
of mid-infrared luminosity profiles. 


We first examine the lack of a significant fraction of counter¬ 
rotating stars in the thick disc. If the thick disc stars were ac¬ 
creted in a succession of minor mergers, we would expect sim¬ 
ilar fractions of prograde and retrograde stars, unlike what we 
observe. If the thick d iscs were created in a few relatively ma¬ 
jor mergers (as in e.g. lRead et al.ll200^ or as a result of the fu - 
sion of a few protogalactic mergers, as in iBrook et"^ (l2004li . 
there would be no guarantee that a particular galaxy would have 
counter-rotating stars. Indeed, in the extreme case where only 
one such an event occurs, there is a 50% chance of having a pro¬ 
grade merger and a 50% chance of having a retrograde merger. 
Thus, the lack of counter-rotating material in ESO-533-4 would 
not rule out a merger origin. If we consid er the ID spectroscopy 
study by lYoachim & Dalcanton ( 2008b|) an d our knowledge of 
the Milky Wav kinematics dchiba & Beeril2000l) . we find that 
none of the five massive galaxies (vc > 120kms“') for which 
thick disc kinematics are known has a counter-rotating disc. 
Therefore, if all thick discs in high-mass galaxies share a sim¬ 
ilar history (which is a reasonable assumption, given that the 
ratio of the masses of the thin and thick discs are roughly sim - 
ilar for most of the high-mass galaxies; IComeron et aLlBOl^ . 
a major merger origin is disfavoured by the currently available 
data. We thus tentatively propose that the thick disc of ESO 533- 
4 has formed either from dynamical heating or that it was born 
dynamically hot in a turbulent disc with a large star formation 
rate. 

To distinguish between those two mechanisms, we can use 
our stellar population data. We see that in most of the tiles in 
Eig. |5] the thin and the thick disc populations are clearly differ¬ 
entiated in the Age - log (Z/Zq) plane, which would go against 


a secular heating scenario where a continuous transition between 
the thin and the thick discs is expected. Instead, it would favour 
a scenario in which the thick disc became thick at high red- 
shift in a relatively fast process. We however caution against 
over-interpreting the stellar population results because of the ex¬ 
tremely dusty conditions under which the fitting is done and be¬ 
cause of their dependence of the smoothing introduced by the 
regularization of the data. Confirmation of our stellar population 
results with longer wavelength spectra is thus desirable. 

Another and probably more robust way to distinguish be¬ 
tween the two mechanisms compatib le with our kinematical data 
is hinted at in iMinchev et al.l (l2015l) . In their model, age gradi¬ 
ents naturally appear in the thick disc-dominated region. These 
gradients would immediately confirm or infirm the formation 
scenario where discs are made at high redshift. Unfortunately, 
because of S /N constrains, we can only make a single bin out 
of the entire thick disc dominated region, which prevents us 
from testing for the presence of an age gradient. Even if a 
larger nu mber of thick d i sc bin s could be made, the effect pre¬ 
dicted by IMinchev et al.l (12015h would be very subtle because 
our observations exte nd only out to about 1.5 disc scale lengths 
(IComeron et al.ll2()T^ . whereas a large age difference (of the or¬ 
der of 2 - 4 Gyr) would in theory be observed only over ranges 
of ~ 3 scale lengths. Larger field of view observations (e.g. with 
MUSE) are thus needed to settle this issue definitely. 

The conclusion is that, if we assume that all thick discs in 
high-mass galaxies share a common origin, observations favour 
a dynamical heating or turbulent disc origin for the thick disc in 
ESO 533-4. A minor merger origin can be discarded, whereas 
a major merger origin is possible if a variety of scenarios are 
able to produce a thick disc in high-mass galaxies. The stellar 
population map of ESO 533-4 is affected by dust absorption and 
should be interpreted with caution. Recovering the stellar popu¬ 
lations from a spectrum is an ill-posed problem which requires 
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the choice of a regularization parameter to smooth the solution. 
Having said that, the stellar population fits seem to indicate that 
this particular galaxy has differentiated thin and thick disc stel¬ 
lar populations, as expected from a short event (either in an early 
turbulent disc or in a major merger). It thus seems to disfavour 
a secular evolution origin. A study of the axial gradients of the 
stellar ages in the thick disc is required to confirm this last asser¬ 
tion. 
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